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Abstract Analytic and numeric approximations are studied in detail for a 
hydrodynamic parameterization of single-particle spectra and two-particle 
correlation functions in high energy hadron-proton and heavy ion reac- 
tions. Two very different sets of model parameters are shown to result 
in similarly shaped correlation functions and single particle spectra in a 
rather large region of the momentum space. However, the absolute nor- 
malization of the single-particle spectra is found to be highly sensitive 
to the choice of the model parameters. For data fitting the analytic for- 
mulas are re-phrased in terms of parameters of direct physical meaning, 
like mean transverse flow. The difference between the analytic and nu- 
meric approximations are determined as an analytic function of source 
parameters. 

1 Introduction 

In 1994-95, a series of papers were written by the Buda-Lund collaboration on the 
study of particle correlations and single-particle spectra for non-relativistic, three- 
dimensionally expanding as well as for relativistic, one-dimensionally expanding 
or three-dimensionally expanding finite systems 0, E|. In these papers, it has 
been emphasized for the first time, that observation of the "true" sizes of particle 
sources is possible only if the single-particle spectra and the two-particle correlation 
functions are simultaneously analyzed. The reason was also given: the HBT radii 
(effective sizes measured by correlation techniques) were found to be dominated by 
the shorter of the geometrical and the thermally induced lenght-scales, while the 
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width of the rapidity distribution or the slope of the transverse mass distribution 
is found to be dominated by the longer of the geometrical and the thermal scales. 
The appearance of the thermal lenght-scales is related to flow and temperature 
gradients, i.e. to the change of the mean momentum of the emitted particles with 
changing the coordinates of the particle emission. Within a thermal radius, these 
changes are not bigger than the width of the local momentum distribution. 

This important effect is presently being re-discovered by various other groups, 
as a consequence of the emerging simultaneous analysis of particle correlations and 
spectra, proposed first in 1994 by the Buda-Lund collaboration. The question 
arizes: Is it possible to uniquely determine the geometrical source radii ("true" 
source sizes) from a simultaneous analysis of particle correlations and spectra? This 
question is basically the same as the analogous question in the momentum space: 
Is it possible to uniquely determine the freeze-out temperature and the transverse 
flow from a simultaneous analysis of particle spectra and correlations? 

We prove by an example that if the absolute normalization of particle spectra 
is not given, then it may be impossible to select from among different minima 
based only on the shape of the single-particle momentum-distribution and on the 
two-particle correlation functions. 

We perform the analysis with the help of an analytically as well as numeri- 
cally well studied model, the hydrodynamical parameterization of the Buda-Lund 
collaboration gj. The domain of applicability of the analytical approximations is 
determined numerically in ref . || for this model. In ref. |r| , this model is shown to 
describe the single-particle spectra and the two-particle correlations at (tt/K) +p 
reactions at CERN SPS simultaneously. 

The same model was tested in refs. [| against the preliminary NA44 data 
on S + Pb reactions at CERN SPS, however, it was shown there that it is difficult 
to find unique, reliable values of the fit parameters. The sources of the difficulty 
are the lack of absolute normalization of spectra and the experimental difficulty of 
proper estimate of systematic errors. We know from earlier fitting of this model 
to spectra without absolute normalization || , that the final results are rather sen- 
sitive to errors and normalizations. In fact, we compare here the two physically 
different minima, found by fitting the Buda-Lund hydro model of ref. || to NA44 
preliminary data on S + Pb central reactions at 200 AGeV. In ref. || results were 
reported on fitting simultaneously the recently obtained absolutely normalized but 
still preliminary particle spectra and final correlations data for the same reaction 
as analysed in ref. ||. 

In the next section the hydrodynamic model is presented, along with a new 
reparameterization of the basic formulas. An approximate analytic solution to this 
model is formulated in a new manner. For a comparision, a numerical approxima- 
tion method is also schemed up. In the subsequent section, radius parameters and 
single particle spectra are calculated using the analytic and numeric methods. The 
transverse mass and the rapidity dependence of the results are shown for a sub- 
stantial range of momentum space. The results are also used to make estimations 
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to the systematic errors introduced by the particular approximations. Finally, we 
summarize and emphasize the importance of the experimental determination of the 
absolutely normalized single particle spectra. 



2 The model and its re-parameterization 

The hydrodynamic model of ref. || is briefly recapitulated below, in a general 
form. The analytic results are then reformulated with new notation. A numerical 
evaluation scheme is also summarized afterwards. 



2.1 The model 

The Buda-Lund model || model makes a difference between the central (core) and 
the outskirts (halo) regions of high energy reactions. The pions that are emitted 
from the core consist of two types: a) They could be emitted directly from the 
hadronization of wounded, string-like nucleons, rescattering with a typical 1 fm/c 
scattering time as they flow outwards, b) Alternatively, they could be produced 
from the decays of short-lived resonances such as p, N*, A or K*, whose decay 
time is also of the order of 1-2 fm/c. This core region of the particle source is 
resolvable by Bose-Einstein correlation measurements. In contrast, the halo region 
consists of decay products of long-lived resonances such as the u>, rj, i]' and Kg, 
whose decay time is greater than 20 fm/c. This halo is not resolvable by Bose- 
Einstein measurements with the present techniques, however, it is affecting the 
Bose-Einstein correlation functions by suppressing their strength. 

In general, the following emission function S c (x,p) applies to a hydrodynam- 
ically evolving core of particle source: 

q / w4 9 d 4 ^(x) Pf , 

= w» 7 To*. — ( ) 

where the subscript c refers to the core, the factor d 4 T l fJ, (x)p fl describes the flux of 
particles through a finite, narrow layer of freeze-out hypersurfaces. The statistics 
is encoded by s, Bose-Einstein statistics corresponds to s = —1, Boltzmann ap- 
proximation to s — while the Fcrmi-Dirac statistics corresponds to s = +1. The 
four- momentum reads as p — p^ = (E p ,p). The four-coordinate vector reads as 
x — x^ — (t,r x ,r y ,r z ). For cylindrically symmetric, three-dimensionally expand- 
ing, finite systems it is assumed that any of these layers can be labelled by a unique 
value of r = yjt 2 — r\, and the random variable r is characterized by a probability 
distribution, such that 

d 4 'S fJ '(x)pf 1 = mt cosh [tj — y] H (r)dr To drjdr x dr y . (2) 
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Here mt = y m 2 + p 2 + p 2 stands for the transverse mass, the rapidity y and 

the space-time rapidity r\ are defined as y — 0.5 log [{E + p z )/(E — p z )\ and r\ — 
0.51og \{t + r z )/(t — r z )] and the duration of particle emission is characterized by 
H(t) oc exp(— (t — to) 2 /(2At 2 )). Here To is the mean emission time, At is the du- 
ration of the emission in (proper) time. The four- velocity and the local temperature 
and density profile of the expanding matter is given by 



l (x) — cosh [77] cosh [n t ], smh[?y t ] — , sinh[^] — , sinh[^] cosh[?7 t ] , (3) 

V n n ) 



sinh[?7 t ] = 6—, r t = Jr% + r%, (4) 
t v 

assuming a linear transverse flow profile. The inverse temperature profile is charac- 
terized by the central value and its variance in transverse and temporal direction, 
and we assume a Gaussian shape of the local density distribution: 

(5) 
(6) 



1 


1 


T(x) 


T~o 


fj,(x) 


Mo 


T(x) 


To 




where fi(x) is the chemical potential and T(x) is the local temperature character- 
izing the particle emission. 

2.2 Core/halo correction 

The effective intercept parameter A» (y, mt) of the Bose-Einstein correlation function 
measures the fraction of pions from the core versus the total number of pions at a 
given value of p, when interpreted in the core/halo picture |ITj, |ll|, |l2| ]. With this 
factor the total invariant spectrum in y rapidity and transverse mass m t follows as 

1 f/2 "< ' f S c ( X , P )d*x. (7) 



dy dm 2 yf\l dy dm\ n 

The momentum dependence of A* parameter is to be measured by the experimen- 
tal collaborations. The experimental determination of \*(p) is very important not 
only because it gives a measure of the contribution of the core to the total number 
of particles at a given momentum, but also as it provides a measure of the mean 
transverse flow and a new signal of partial {7a(1) symmetry restoration H. There- 
fore it is strongly recommended that experiments report this A*(p) parameter of 
the Bose-Einstein correlation function and not just present partial fit results, like 
the momentum dependence of the radius parameters. 
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2.3 Re-parameterization 

The original version of the core model contains 3 dimensionless parameters, a, b and 
d, that control the transverse decrease of the temperature field, the strength of the 
(linear) transverse flow profile and the temporal changes of the temperature field, 
respectively, keeping only the mean and the variances of the inverse temperature 
distributions. These are very useful in obtaining simple formulas, however, they 
make the interpretation of the fit results less transparent. Hence we re-express 
them with new parameters with more direct physical meaning. 

The surface temperature is introduced as T r — T(r x — r y — Rq, t — r ) and 
the "post-freeze-out" temperature denotes the local temperature after most of the 
freeze-out process is over, T t = T(r x = r y = 0;r = r + \/2At). Here Rq stands 
for the transverse geometrical radius of the source, tq denotes the mean freeze-out 
time, At is the duration of the particle emission and we denote the temperature 
field by T(x). The central temperature at mean freeze-out time is denoted by 
T = T(r x = r y = 0; r = to). 

Then the relative transverse and temporal temperature decrease can be intro- 
duced as 

(AT/T) r = (8) 

1 r 

(AT/T) t = ^-Z*, (9) 
J-t 

and it is worthwhile to introduce the mean transverse flow as the transverse flow 
at the geometrical radius as 

(u t ) = b^. (10) 

TO 

The dimensionless model parameters can thus be expressed with these new, physi- 
cally more straightforward parameters as 



a 2 



R G 



7T (AT/2%, (11) 



b = (u t ), (12) 



Rg 
At 2 



(AT/2%. (13) 



Note, that eqs. (||) and ([lO]) were introduced earlier in ref. 0] also to simplify 
the interpretation of data fitting. We present the complete re-parameterization 
herewith, including egs. ([)[) and the re-parameterization of both the radius param- 
eters and single-particle spectra. 
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2.4 Analytic approximations 

In Ref . H , the Boltzmann approximation to the above emission function was eval- 
uated in an analytical manner, applying approximations around the saddle point 
of the emission function. The resulting formulas express the Invariant Momentum 
Distribution (IMD) and the Bose-Einstein correlation function (BECF). Now we 
re-express the formulas given in ref. j| with the help of our new parameters. 
The particle spectra can be expressed in the following simple form: 



jV(p) = * £FCexp(- £ )> ( 14 ) 



(2tt) 3 P ^ T{x) 

E = m t cosh(?7), (15) 

V = ^L, (16) 



A 



T 



Here the quantity x stands for the average value of the space-time four-vector pa- 
rameterized by (r, fj, r^,r^), denoting longitudinal proper-time, space-time rapidity 
and transverse directions. These values are given as 

T = T , (18) 

77 = ; A 2 "» t ' (19) 

Jo 

= (u t )R G = — , (20) 

rv =0. (21) 

E denotes the energy of a particle from the center of particle emission, measured 
in the Longitudinal Center of Mass System (LCMS) frame. The effective volume is 
denoted by V, see below for details, and the correction factor C takes into account 
the effects of long-lived resonances and the deviation of the saddle-point result 
from the more possible naive expectation, which would be the same expression 
with (7=1. The notation a denotes an invariant quantity a, that depends on 
y — j/o j m t, T and the other parameters of the model in a boost-invariant manner. 
In ref. |j| this was denoted by a^, for example, Arj in the present paper was de- 
noted by At7*. The average invariant volume of particle emission in eqs. ( |l4| , |l6| ) is 
given by V that is in our case the product of the average (momentum-dependent) 
transverse area (27ri?j_ ), the average (momentum-dependent) longitudinal source 
size (27r) 1 / 2 j R || . The averaging of the size of the effective volume over the duration 
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At 

of particle emission is expressed by the factor — ; — . These quantities read as 

At 



At 



Rw = 



At 2 



1 + (AT/T) r 
An 2 



E 



l + An 2 



An 



E 



R 2 



G 



+ ((u t ) 2 + (AT/T) r ) 



E 
To 



(22) 

(23) 

(24) 
(25) 

(26) 



This completes the specification of the shape of particle spectrum. These results 
correspond to the equations given in ref. |jj although they are re-expressed with 
new combination of the variables. Please note that the Boltzmann- factor can be 
expressed approximately as 



exp 



p ■ u(x) — fi(x) 



exp 



Mo 



(y - yo) 2 



exp 



To 



x exp 



T 2{An 2 + To/m t ) 

{u t ) 2 {m 2 t -m 2 ) 
2T [T + m t ((u t ) 2 + (AT/T) r )] 



(27) 



The HBT radius parameters were evaluated in ref. || in the following way: the 
space-time rapidity rj s was defined as the solution of the equation 



dS(n s ) 
drj 



(28) 



The Longitudinal Saddle Point System (LSPS) was introduced as the frame where 
rj s = 0. At the so-called saddle-point one has 



dS_ 
~dr~ 



dS 
drj 



dS 

dr x 



as 

dr. 



= 0. 



(29) 



This resulted in a set of transcendental equations for the position of the saddle- 
point. These equations for the positions were solved approximatelly with the help 
of an expansion of the transcendental equation in terms of small parameters, like 
the deviation of rj s from the mean rapidity of the pair. The calculation resulted in 
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the following value for rj s : 

yo-y 

Vs = y+ mt :■ ( 3 °) 

l + Ar^-* -1) 
Jo 

Hence the relative rapidity of the LSPS frame as compared to the LCMS frame is 

l yo-y , Q1 s 

l + A^-^-l) 
Jo 

Note that LCMS is the frame where the rapidity belonging to the mean momentum 
of the pair vanishes, y — ]l3]j . In ref. S this quantity was denoted by the slightly 
more complicated notation 

^CMS = V L } 

and the maximum of the Boltzmann 
factor, which in our present notation reads as rj = ryj and stands for a modified 
saddle-point, used for the calculation of the particle spectra in ref. ||. Note that 
may deviate from rj substantially at low values of m t , especially if T > m. 
This frame, defined by the maximum of the Boltzmann factor, given by rj 
in the LCMS, plays a key role in the calculations. Hence this frame, introduced 
in ref. [[| without a name, deserves a name. We suggest Longitudinal Boltzmann 
Center System (LBCS). In general, this LBCS frame is defined by the solution of 
the 

= (32) 

on 

equation, where /b = exp (— [pu(x) — [i(x)]/T(x)) stands for the Boltzmann factor 
only, but does not include the Cooper-Frye flux term. A more detailed comparision 
between the LBCS and LSPS frames is in preparation for a separate publication jl4j ]. 
In ref. || the spectrum was evaluated in the LBCS frame, while the correlation 
functions in the LSPS frame. Let us recapitulate the results for the two-particle 
correlation functions: The mean momentum is denoted by K = (pi +p2)/2, its 
components are (Kq, Kl, K t , 0) in the (i, r z ,r x ,r y ) reference frame. The transverse 
velocity (3 t reads as 

3 _K t _K t 1 _ V t 

' K M t cosh(y) cosh(y)' v ; 



The effective source parameters are obtained as 

1 1 1 



(34) 



A *l* A?7 {fj -> v l ) cosh 2 (t]L ) ' 

Rj = r^nl = r s 2 A^ 2 (7J^^), (35) 

R\ = Rl = Rl 2 (r}^^), (36) 

R 2 T = Ar 2 = Ar 2 (rj -> 77^). (37) 
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The modified position of the maximal emissivity is given in the LCMS frame by 
(T s ,rj s ,r XjS ,0), where 

r x , s = r^(r]->Vs), (38) 

T s = T = TO. (39) 

This notation means that, for example, Ary, corresponds to the function Arj when 
the variable rj is replaced by the quantity rj^. Note, that A77* corresponds to the 
effective space-time rapidity width of the particle emission function in an LSPS 
calculation, while A?y corresponds to the width of the Boltzmann-factor only in the 
LBCS frame. Similar relations hold for R and At and At*. The Bose-Einstein 
correlation functions can be written in the so-called Bertsch-Pratt side-out-long ref. 
frame as 



C(Ak,K) = 



r>2 
Ft 2 



Rl 



p2 



1 + A* CXp ( — R 2 s ideQ 2 s ide ~ RoutQout ~ R\Q"l — Rout,hQ hQout) j 

(40) 

R\ = Rl (41) 

(42) 



E»2 _r_ C d2 



cosh (y) 
1 

cosh 2 (y) 
Vt 

cosh 2 (y) 



cosh 2 (?7 s )i? 2 + sinh 2 (?7 s )i? ; 



cosh 2 (77^)i? 2 + sinh 2 ^).?? 2 . 



cosh(ry s ) sinh(?7f )i? 2 + sinh^) cosh(ry^)i? 



(43) 
(44) 
(45) 



Here we have utilized the form of equations given in ref. [^| and the transformation 
of K L to M t cosh(y) as introduced in ref. . 

Although the two-particle Bose-Einstein correlation function is manifestly co- 
variant [ pl| , its Bertsch-Pratt parameterization is frame-dependent. A possible 
covariant parameterization ]17| , |lSj , applied recently to the particle interferometry 
by the Regensburg group |19|, [2C|. We find that the simplest possible covariant 
generalization is not exactly the YKP parameterization, but the formulation given 
by the Buda-Lund collaboration in ref. ||, see especially eqs. (44) and (21-26) of 
ref. H . Here we repeat only the results after Gaussian approximation to the source 
function: 



C{Ak,K) = 1 + A» exp(-Q 2 i? 2 - Q 2 i?f - Q t 2 i? 2 ), (46) 

Qt = QLe + QL: (47) 

Q T = Q ■ n{x L s ) = Q cosh(r^) - Q z sinh(^), (48) 

Qn = yfQ-Q-(Q- n( X t)Y -Q 2 t =Qo sinh(r/ s L ) - Q z coshfaf). (49) 
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where n(x^) — (cosh(ryf'), 0, 0, sinh(ryf')) is a normal-vector at a; s in LCMS. Note, 
that this formulation is equivalent with the YKP formulation, however the correla- 
tion function is given by a simple purely quadratic form in the present formulation, 
in contrast to the YKP expression, where the invariant relative momentum com- 
bination Q v is written out explicitely in terms of its non-invariant components. A 
more detailed comparision between the YKP and the Buda-Lund parameterization 
is being discussed in ref. fl4|| . 

2.5 Numeric approximations 

Approximate single particle spectra and Bose-Einstein correlation functions can be 
calculated by numerical integration of equation (1), as well, expressing the means 
and the variances of the hydrodynamically evolving core of particle emission. The 
original method of 'means and variances' was proposed first by the Regensburg 
group in ref. [lq ]. We shall utilize herewith the core/halo corrected version of 
these relations as given recently in ref. [ pT| . The limitation of these approximations 
is discussed in refs. || . For example, possible double-Gaussian structures or non- 
Gaussian features are neglected in this approximation. The analytic approximation 
yields Gaussian functions in proper time r and space-time rapidity 77, hence includes 
a deviation from Gaussian shape in t and z. The numerical approximation assumes 
Gaussian forms in t and z that is not so well suited to the kinematics of ultra- 
relativistic reactions as Gaussians in r and n. The spectra and the HBT radius 
parameters are defined in the Gaussian core/halo model approximation as follows: 



C(K,Ak) = 1 + A*(K) exp(-.R? j (K)Ak i Ak i ) ) (50) 

A*(K) = [Y C (K)/A(K)] 2 , (51) 

i^.(K) = {{x i -p i t){x j -p j t)) c -{{x i -p i t)) c {{x j -p j t)) c , (52) 

(f(x,K)) c = / dV(a:,K)S c (a:,K), (53) 



where i,j — side, out or long as before, and S c (x, K) is the emission function that 
characterizes the central core, as given by eq. The mean momentum is defined as 
K = 0.5(pi + P2), the relative momentum is given by Ak = pi — P2. The spectra 
of all the bosons and the spectra of the bosons from the core described as 

N(p) = (1) = J d 4 x[S c (x 1 p) + S halo (x lP )] 1 (54) 
N(p) c = (l) c = J d 4 xS c {x,p). (55) 

In this picture, the reduction of the intercept parameter is the only effect on the 
correlation function that stems from the halo, the variances of the core correspond to 
the Gaussian core/halo model radii of the measured correlation function. Although 



On source parameters 



11 



the above expressions are formally similar to the original version of Gaussian model- 
independent radii of ref. |lq |, they cannot be obtained with an expansion around 
Q = 0, as they correspond to a large Q expansion of the Bosc-Einstein correlation 
function |§ 

3 Calculating observables 

Both analytic and numeric approximations were used to calculate and show the mo- 
mentum dependence of the observables from the hydrodynamical model of ref. || . 
These observables are the effective radius parameters (HBT radii), the shape and 
the slope of the single-particle spectra. Figures are drawn for two very different 
sets of source parameters (or model parameters). To recall, the particle source is 
characterized by the means and the variances of the density distribution, the inverse 
temperature distribution and a linear flow. This yields 9 free parameters, /xo,To, 
r , R G , At, A V , (AT/T) r , (u t ), (AT/T) t , respectively. 

The examined momentum space is divided into 40x40 sub-intervals in mt and 
y dimensions, respectively, that allowed for fine resolution of the distributions. As 
a drawback, with such a resolution the numeric integration version takes much 
longer time than the analytic one. In the presented case generation of data took 
10 hours for one run. Scales on the pictures are kept the same for the same kind of 
distributions except for the average emisson rate that differs remarkably for the two 
sets of model parameters. The actual parameter set values (with = 0) used in 
the particular calculations are indicated below each drawing and they are denoted 
by names Source Parameters Set 1 and Source Parameter Set 2. Parameter Set 
1 was obtained from a ref. (^] , fitting absolutely normalized spectra in the NA44 
acceptance, while Parameter Set 2 was obtanined in ref. fitting unnormalized 
preliminary particle spectra together with correlation data. Note that one can 
distinguish between different model parameters only if the value of no is known 
from other observations than the ones already exploited in the present analysis. 

3.1 Analytic results 

The two model parameter sets mentioned above were applied to the analytic ex- 
pressions formulated in the previous section by eqs. (|l4|) to ([i"5|). See Figures 1 
and 2 for details of the momentum space distributions of the observables. Notice 
the substantial difference of the particle spectra for Parameter Set 1 and Parameter 
Set 2. Also notice the deviations of the radius parameters at small mt and at large 
relative rapidities to midrapidity for the two different source parameter sets. Along 
with the comparision to the numeric results later this reflects the limitations of this 
kind of approximation. See ref. |14j for an improved treatement. 
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3.2 Numeric results 

The numerically evaluated HBT radius parameters and single particle spectra are 
obtained from the eqs. (^0U52|j5^ , |55| ) , utilizing the Boltzmann approximation to the 
source function of eq. (^J). Note that this scheme is not an exact calculation, but an 
approximation in a different way than the analytic approach, therefore it is suitable 
to estimate the systematic errors of the model parameters and to cross-check the 
uniqueness of the minimum in fitting the model to experimental data. See Figures 
3 and 4 for details of HBT radius parameter distributions as well as single particle 
spectra in this approximation scheme. 

3.3 Differencies between the analytic and the numeric results 

The differencies between the observables as calculated from the two sorts of model 
approximations are presented on Figures 5 and 6. From these drawings one can 
learn the critical ranges where the two approximation schemes differ from each 
other beyond a given tolerance level. 

Figures 5 and 6 show the rapidity and transverse mass dependence of the rela- 
tive deviations between the analytic and the numeric approximations. For example, 
let us consider the top left panel of Figure 5. On this panel, the relative deviation 
of the numerical and analytic approximations for the side radius component is eval- 
uated in the following manner: The analytical result for a (y,m t ) bin i is denoted 
by di, the numerical result is denoted by rij. Then the relative deviation between 
the two approximation schemes is defined as 

5 2 a t {ai - rii) 2 

— = 2 > ( 56 ) 

which is plotted on the top left panel for the side radius parameter and in subsequent 
panels for the out, long and cross term, the spectra and the slope of the spectra on 
the subsequent panels. This quantity will be understood and estimated in the next 
subsection as a function of some small expansion parameters, that are analytically 
obtainable for any set of source parameters. In turn, this result can be utilized to 
improve the fits of the analytic expressions to measured data. 

3.4 Estimating systematic errors of approximations 

The aim of the present subsection is to analytically understand the systematic 
errors on the analytic approximations that we utilize to evaluate the spectra and 
the HBT radius parameters. 

In a fit with the analytic approximations, the % 2 of the fit is given as 
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where di denotes the measured data point at a given bin, for example, R s ide(yi> m t,i), 
and the experimental error on this quantity is given by e%. A numeric fit to the 
same data minimizes the following 

numeric \n distribution. This may have different minima for the model parame- 
ters from the minima of \\. However, the x n distribution can be approximately 
reconstructed from an analytic fit, as follows: 

Xl - Xt + Sxl, (59) 

Sxl = £^=£^4, (60) 

% 1 % * 

where 5 2 a.i — (ctj — rii) 2 is the difference between the analytic and numeric result, 
and can be regarded as the systematic error of the analytic approximation, while 
the relative systematic error of the analytic calculation is S 2 Ai = ^-jjK (Keep in 

a i 

mind that dj can be any of the analytically evaluated radius parameters or analytic 
result for the particle spectra). Our purpose is to obtain approximate analytical 
expressions for the relative error of the analytical approximations, <5 2 A;. These 
quantities are shown in Figures 5 and 6. 

Figures 5 and 6 indicate large relative errors in certain regions of the (y,m t ) 
plane. These regions coincide with the regions where the so called small expansion 
parameters || of the model start to reach values close to 1. The analytic expres- 
sions for the observables (HBT radius parameters and single-particle spectra) were 
obtained in ref. Q under the condition that the parameters (77^, A77*, r^^/ro) are 
all much less than 1. This was due to the approximate nature of the solution of 
the saddle-point equations, and the expansion of the transcendental equations in 
terms of these small parameters. Figures 7 and 8 show the (y, m t ) dependence of 
these small parameters and their squares. We observe the expected similarities to 
the distributions of the relative errors 6 2 Ai, Figure 5 and 6. The different small 
parameters become large in well separated domains of the momentum space (typ- 
ically below lOOMeV and above 1 GeV), e.g. 77^ at small mt and large \y — yo\, 
A77* at small mt and small \y — j/o I ; r XtS /TQ at large m t independently of y. As a 
consequence, the relative errors of the analytic approximations for any radius pa- 
rameter or momentum distribution can be parameterized as a linear combination 
of the squared small parameters: 

S 2 A = c?( V f) 2 + + ci(r x , s /r Q ) 2 . (61) 

The three constants {cf , c^,c^) are determined for the observables A = R S idei Rout, 
Riong, Ro U t-iong> -^(p)j T e f f{y,m t ) . Since the small parameters are expressed as 
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a function of the model parameters or "true" source parameters, eqs. (p5|-|37|) , the 
desired analytical formula for all the systematic errors of the evaluation of all these 
6 observables is given in the form of eq. (|6l]) . 

Note that the region, where the analytic expressions are most precise, corre- 
sponds to a curved region in the (y, mt) plane, that at low m* starts off-mid-rapidity 
and at high shifts to mid-rapidity, in case of pions. This region almost exactly 
coincides with NA44 acceptance for pions. For heavier particles all the 3 small 
parameters decrease substantially. The analytic calculation is thus more precise 
for heavier particles than for pions. 

For determining the numerical coefficients cf, the distributions of the differ- 
encies on Figures 5 and 6 the CERN optimizing package MINUIT was used 
MINUIT finds the minimum value of a multi-parameter function and it analyzes 
the shape of the function around the minimum. The fitted error distributions in 
terms of the small parameters are shown on Figures 9 to 12 together with the best 
estimates of the systematic errors with the help of eq. [5l[ The coefficients cf of 
the parameterized systematic error distributions along with their errors are shown 
in Table 1 and Table 2. 

Note that the coefficients cf are found to be only weakly dependent on the 
source parameter sets for most of the cases, they are all smaller then the coefficents 
for Parameter Set 1 multiplied by 2. On Figure 11, the fit of (S N (p) / N (p)) 2 reflects 
the effect that above 0.5 MeV the calculation was forced to take the amplitudes 
with low weight due to the fact that the absolute values of N(p) are very small in 
this range. 



4 Conclusion 

The first combined HBT and spectrum analysis, reported by the Buda-Lund col- 
laboration in 1994-95, is re- visited here for the purpose of a systematic numerical 
and analytical evaluation of the model. We identified the regions in the rapidity - 
transverse mass plane where the analytic approximations and the numerical ones 
deviate from each other. These regions were found to coincide with the regions 
where the small expansion parameters of the analytic approximation start to grow 
significantly. The deviation between the analytical and the numerical results is 
characterized by positive definite quadratic polinomials built up from the small 
parameters. We find that the NA44 acceptance is ideally suited for the precise 
evaluation of the Buda-Lund model. 

As a by-product we find that the parameters of the Bose-Einstein correlation 
functions as well as the shape and the slope parameters of the double differential 
invariant momentum distribution are similar within 10% for two physically very 
different model parameter value sets in the momentum space domain where the 
approximations are valid. However, these sets result in a factor of 7 - 10 difference 
in the absolute normalization of the single-particle spectra. Hence, it is strongly 
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recommended to publish the experimentally measured single-particle spectra with 
their absolute normalization for future CERN and RHIC heavy ion experiments, 
for as many type of particles as possible, as a two-dimentional function of y,mt- 
Figures 1 to 4 illustrate that the y and the mj dependence of N(y,m t ) can not be 
factorized. 

Of course, our example of two physically different parameter sets resulting in 
similar correlations and unnormalized spectra does not imply that such similarity 
is achieved for any two different parameter sets. The role of each parameter can 
be investigated, for instance, analytically like in ref. j| ||. The typical behaviour 
that one expects is that for different values of the model parameters the particle 
correlations and spectra are different. We would like to remind the readers that 
the existence of certain scaling limiting cases was pointed out also in ref. || , where 
the dependence of the HBT radii on some of the model parameters was analytically 
shown to vanish in certain domains of the parameter space. 

Note that the intercept parameter of the correlation function is included in 
the core/halo correction factor, hence it is also strongly recommended to publish 
the experimental HBT results including a momentum-dependent determination of 
A*, too. 
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Table captions 

Table 1. Coefficients cf of the parameterized systematic error distributions for Source 
Parameter Set 1. 

Table 2. Coefficients cf of the parameterized systematic error distributions for Source 
Parameter Set 2. 

Figure captions 

Fig. 1. Simultaneous results for particle spectra and HBT radius parameters. The 
analytic approximations were utilized to evaluate the model for Parameter 
Set 1. 

Fig. 2. Simultaneous results for particle spectra and HBT radius parameters. The 
analytic approximations were utilized to evaluate the model for Parameter 
Set 2. 

Fig. 3. Simultaneous results for particle spectra and HBT radius parameters. The 
numeric approximations were utilized to evaluate the model for Parameter 
Set 1. 

Fig. 4. Simultaneous results for particle spectra and HBT radius parameters. The 
numeric approximations were utilized to evaluate the model for Parameter 
Set 2. 

Fig. 5. Relative deviations between the analytic and the numeric approximations to 
the HBT parameters and particle spectra for Parameter Set 1. 

Fig. 6. Relative deviations between the analytic and the numeric approximations to 
the HBT parameters and particle spectra for Parameter Set 2. 

Fig. 7. Small expansion parameters of the hydrodynamic core model for Parameter 
Set 1. 

Fig. 8. Small expansion parameters of the hydrodynamic core model for Parameter 
Set 2. 

Fig. 9. Paramctrization of the relative deviations between the analytic and the nu- 
meric approximations to the HBT parameters and particle spectra for Pa- 
rameter Set 1. 

Fig. 10. Parametrization of the relative deviations between the analytic and the nu- 
meric approximations to the HBT parameters and particle spectra for Pa- 
rameter Set 2. 

Fig. 11. Parametrization of the relative deviations between the analytic and the nu- 
meric approximations to the HBT parameters and particle spectra for Pa- 
rameter Set 1. 
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Fig. 12. Parametrization of the relative deviations between the analytic and the nu- 
meric approximations to the HBT parameters and particle spectra for Pa- 
rameter Set 2. 
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S 2 A ~ 

(5R s ide/Rsidej 2 
(6 R out /R out ) 
(SRlong/ Rlong) 
( SR lutl/ R lutl) 2 

(SN(p)/N(p)f 
(ST eff /T eff f 



0.015 ± 0.001 
0.013 ± 0.001 
0.027 ± 0.001 
0.130 ± 0.001 
0.014 ± 0.001 
0.006 ± 0.001 



0.003 ± 0.001 
0.001 ± 0.001 
0.014 ± 0.001 
0.208 ± 0.001 
0.044 ± 0.001 
0.042 ± 0.001 



0.002 ± 0.001 
0.039 ± 0.001 
0.032 ± 0.001 
0.001 ± 0.001 
0.015 ± 0.001 
0.001 ± 0.001 



Tabic 1: 



6*A 



(SR slde /R slde ) 2 0.001 ± 0.001 0.007 ± 0.001 0.003 ± 0.001 

(SRout/Rout) 2 0.008 ± 0.001 0.002 ± 0.001 0.001 ± 0.001 

(SRlong /Rlong) 2 0.018 ± 0.001 0.024 ± 0.001 0.031 ± 0.001 

iSR 2 ouU /R 2 outl ) 2 0.154 ± 0.003 0.398 ± 0.002 0.001 ± 0.001 

l5N(p)/N(p)) 2 0.035 ± 0.001 0.052 ± 0.001 0.001 ± 0.001 

(ST ef f /T ef f ) 2 0.012 ± 0.001 0.042 ± 0.001 0.001 ± 0.001 



Table 2: 
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Fig. 1. 
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Fig. 3. 
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Fig. 5. 
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Fig. 7. 
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Fig. 9. 
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Fig. 11. 



source parameters 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9810341vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9810341vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9810341vl 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9810341vl 



This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9810341vl 



This figure "figl-6.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9810341vl 



This figure "figl-7.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9810341vl 



This figure "figl-8.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9810341vl 



This figure "figl-9.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9810341vl 



This figure "figl-lO.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9810341vl 



This figure "figl-ll.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9810341vl 



This figure "figl-12.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9810341vl 



